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Abstract 

We study electronic properties of solids with correlated d electrons which could 
be described by a multiband Hubbard Hamiltonian in the weak-interaction case, 
U/w < 1. The one-electron part of the many-body Hamiltonian is described by a 
tight-binding LMTO method. The many-body part is treated by non-selfconsistent 
FLEX- type approximations with adjusted chemical potential not to change the 
LMTO band filling. The calculated DOS gets narrower but is only little influenced 
by electron correlations at the Fermi energy. A precursor of a satellite band in the 
paramagnetic bcc Ni is found at about 6 eV below the Fermi level in agreement 
with experiment. 



INTRODUCTION 

A general scheme for determination of the electronic structure of solids with corre- 
lated electrons is still missing. The density functional formalism within the local density 
approximation (LDA), in which the electron-correlation part is guessed on the basis of 
the homogeneous electron gas model, has proved to be a highly reliable method for the 
evaluation of the ground state properties On the other hand, there are many exam- 
ples of the failure of the LDA, e.g., the study of the excitation spectra of solids or the 
evaluation of the gap in insulators and semiconductors. Another example of the failure of 
the LDA are solids whose electronic structure is better described in terms of atomic-like 
electronic states rather than in terms of the homogeneous electron gas model on which 
the LDA is based. 

The electronic properties of strongly correlated materials which cannot be ade- 
quately described within the LDA are usually studied in the framework of simplified 
models like, e.g., the single-band Hubbard model (for a recent review see Ref.0). The 
correlation effects in solids are usually classified in terms of the ratio of the on-site Coulomb 
energy U to the bandwidth w. One distiguishes three regimes, namely, (i) the weak inter- 
action case (U/w < 1, transition metals), (ii) the intermediate interaction case (U/w ~ 1, 
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metal-insulator transition regime and Kondo systems), and (iii) the strong interaction 
case (U/w > 1, rare-earth systems and wide gap solids). It should be mentioned that a 
general approach bridging all the above mentioned cases is still missing even at the model 
level. 

Recently some progress has been made in the reformulation of the LDA (the time- 
dependent LDA j3|) which is useful in the study of excitation energies. The so-called 
GW approximation (GWA) [[| was successfully applied to the gap problem in a number 
of semiconductors and insulators. An alternative approach to the gap problem which is 
based on a straightforward generalization of the LDA is the so-called LDA+U method 
0. The GWA was also applied to metallic systems. Although computationally highly 
demanding, it yielded band narrowing in the photoemission spectra of nickel |J, but it 
failed to produce the well-known satellite structure below the main peak The presence 
of the satellite in Ni was succesfully explained in the framework of the T-matrix formalism 
developed recently on ab initio level (see Ref. || and references therein for previous 
semiempirical approaches to this problem). An alternative approach to the satellite and 
the band narrowing in nickel is based on a three-body scattering approximation that 
employs the Faddeev equations || [10| . 

In a recent paper [llj the so-called LDA++ method was introduced which extends 
the LDA+U method by accounting for dynamical electron correlations. The starting 
point of the LDA++ method is a multiband Hubbard Hamiltonian whose one-electron 
part is identified with the LDA Hamiltonian with a subtracted double counting correction 
for the average Coulomb interactions in the LDA. The LDA++ approximates solutions in 
different correlation limits by different many-body approximations, namely it employs the 
so-called Hubbard I solution fT2| in the strong- interaction limit, the iterated perturbation 
theory (IPT) within the dynamical mean-field theory (DMFT) || for the intermediate- 
interaction case, and the so-called fluctuation exchange approximation (FLEX) [TJJ in the 
weak- interaction limit. 

In the present paper we follow the basic approach of the LDA++ but with some 
conceptual as well as computational differences. In particular, the main effort is directed 
towards a formulation of the problem in a way suitable for future generalizations to ran- 
dom alloys and their surfaces and interfaces. This means that the formalism is strictly 
based on Green functions which allow for configurational averaging in the case of alloys. 
The starting point is a multiband Hubbard Hamiltonian (MBHH) whose parameters are 
naturally determined from the corresponding tight-binding linear muffin-tin orbital (TB- 
LMTO) Hamiltonian (for a review of the TB-LMTO method and its applications to 
random surfaces and interfaces see a recent book, Ref. [14}j). The basic approximation 
adopted here is the assumption on the local (site-diagonal, or wave-vector independent) 
selfenergy which is reasonably well justified for transition metals and their alloys JT5] . 
The local approximation for the selfenergy is also required from the formal point of view 
for the extension of the theory to the case of random alloys within the coherent potential 
approximation (CPA) fl6|| . We employ the single-channel approximations (FLEX) of the 
canonical perturbation theory to solve the many-body part of the problem. 



THEORY 

The electronic structure determined within the LDA is described by a TB-LMTO 
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Hamiltonian with the overlap matrix equal to unity 



^lmto = E |RA>i^ R °, (R'A'I, (l) 

RA.R'A' 

where 

fjLMTO _ r r , a 1/2 q7 a 1/2 / Q x 

^RA,R'A' - ^RAA' "RR' + a ra ^RA.R'A' A R'A' ■ \ Z ) 

Here, R is the site index, A = {La) = {ima) is the spinorbital index, L = (£m) is the 
orbital index, a is the z-projection of the spin, C, A, and 7 are site-diagonal matrices of 
potential parameters and S" 7 is the matrix of structure constants in the orthogonal LMTO 
representation 



•Sra.ra' - - 1 ]ra,r'A' 



(3) 

RA,R'A' 



Here S° is a matrix of canonical structure constants, S 13 is a matrix of screened structure 
constants, and (3 is a site-diagonal matrix of screening constants. 

The parameters of a multiband Hubbard Hamiltonian in second quantization with 
creation (a RA ) and destruction a R A operators 

# HUbb = *RA,R'A' ARA OR'A' + E ^W^RA^RA' (4) 

RA,R'A' R,A,A' 

are found from the LDA calculations. The density operator n R A = a RA a R A- The hopping 
integrals 

*RA,R'A' = ^RA SrA,R'A' ^R'A' > ( 5 ) 

are identified with the site off-diagonal elements of ([!]), while the atomic levels £ra,ra' = 
£raa' = e RA^A,A' are determined from the condition that the average occupancy n R A of the 
state |RA) as calculated within the LDA and within the Hartree-Fock approximation for 
multiband Hubbard Hamiltonian M) are identical. It means 



£ra 



Cra — ^ RA A — Cra — ^2 Urw ^ra' • (6) 



The Hubbard interaction parameter £7 R aa' is determined within the LDA usually as a 
second derivative of E tot with respect to occupancies of the interacting states A, A'. Due 
to the Pauli principle, £/raa' = C^raa'(1 — <5aa')- Here we consider translationally invariant 
solids with one atom per unit cell. The hopping integrals then depend only on the differ- 
ence R — R'. The on-site levels €ra as well as the pair interactions C/raa' are independent 
of R. 

We study four different approximations of the canonical many-body perturbation 
theory, namely, (i) the second order perturbation theory (SOPT), and three approxi- 
mations of infinite order that correspond to electron-electron (hole) scatterings in a sin- 
gle channel, namely, (ii) the T-matrix approximation (TMA) which describes electron- 
electron scatterings, (iii) the random-phase approximation (RPA) corresponding to electron- 
hole scatterings, and (iv) the GW-approximation (GWA) in which the bare interaction 
U is replaced by a (complex and energy-dependent) interaction W{E) (double wavy line) 
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renormalized by repeated excitations of electron- hole pairs (polarization bubbles). The 
graphs for one-particle selfenergy are shown in Fig. 1. We assume 
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2nd order 




TMA 






T,i, 



GWA 




Figure 1. Graphs for one-particle selfenergy in second order perturbation theory (SOPT), T-matrix 
approximation (TMA), random-phase approximation (RPA), and GW-approximation (GWA). 

a local approximation, i.e., the selfenergy is independent of the k- vector, which means 
that it is diagonal in the site representation, Xra,r'a' = ^r,aa'^r,r'- We use the technique 
of causal Green functions at T = in our calculations. It follows from the adopted 
construction of the on-site levels in (^j) that the Hartree-Fock part of the one-electron 
selfenergy is already contained in the atomic levels €ra- Consequently, the expansions of 
the selfenergy start from the terms quadratic in the Hubbard U. The unperturbed one- 
particle Green function corresponding to the Hartree-Fock approximation (HFA) applied 
to H Hubh is identical with the Green function of H LMT °. 

We assume a cubic symmetry of the lattice and confine ourselves to s, p, and d 
states. In addition, we neglect the pair correlations in s and p states and for simplicity 
we assume U\y — U(l — £aa')<^,2, i-e., we neglect exchange interactions and consider only 
the pair interactions between the d-states. The Green functions and selfenergy are then 
diagonal in the spinorbital indices A. The selfenergy can acquire only two values that 
correspond to e g and t 2g representations of the cubic point group. 

As an example we give here the derivation of basic equations for the TMA which 
is an extension of the derivation given in Ref. |17] to the multiband case. The basic 
equations of the TMA in the language of causal quantities (superscript c) read 



duo 
2^ 



Gf{E-u)Gf!{u>) 



U 



U . 



(7) 
(8) 
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-*4 C) (£) = E("1) T ^-TH(E + u)G^(u), (9) 

x , J-oo Z7T 

where ^ x c l,(E) is the two-particle propagator, t[ c x \(E) is the electron-electron T- matrix, 
and Z x c) (E) is the one-electron selfenergy. The factor (—1) in @ is due to a closed fermion 
loop. In numerical calculations, it is advantageous to replace the causal quantities by the 
retarded ones (without superscript) as follows 

= Re#AA'(£) + *sgn(£ - 2/i) Im* AA ,(£) , (10) 

T$(E) = ReT xx ,(E) + % sgn(£ - 2/i) ImT AA ,(£) , (11) 
E x c) (E) = Re£ A (£) + i sgn(£ - /x) Im£ A (£) , (12) 
G A c) (£) = ReG A (£) + isgn(E - fi) lmG x (E) . (13) 

By inserting (|T0l)-(|l"3|) into (0)-(0) we find equations for the retarded quantities which are 
holomorphic in the upper halfplane of the complex energy z and vanish for z — > 0. It is 
therefore sufficient to calculate only their imaginary parts because the real parts can be 
found using a dispersion relation 

,~ da; 

J-oo 7T 0> — 2 

which allows to determine the real and imaginary parts of X(z) for complex z as well as 
for real values of the energy. The equations for the retarded quantities read 



Im* 



AA' 



da; 

(E) = -sgn(£ - 2/i) / — lmG\(E - oj) lmGy(cu) , (15) 

JE-n 71 



Im£ A (£) = - E r — ImT AA'(^ + w) ImGv(^) • (17) 

y J2fi-E 71 



T ^ = i-u*ME) - U ' (16) 

>2fi-E 71 

In the derivation of ( |T5f )- (|l7|) we have used analytic properties of the retarded quantities 
from which follow useful relations between convolutions of real and imaginary parts, as, 
for example 

[°° ^ ReG\(E - to) ReG* M = - [°° ^ \m.G x {E - to) lmG y {uj) . (18) 

The equations for the other approximations can be derived in a similar way. An 
expression for the renormalized pair interaction W(E) requires matrix inversion even in 
the simplified case U X y —U(l — 5 X y) assumed here, since 

E i 6 ^" " U X y,®y,(E)} Wy>y{E) = U X y , (19) 
A" 

where $>\(E) is the electron- hole bubble and 

f da; 

lm^ x (E) = - — lmG x (E + a;) lmG x (uj) . (20) 



The retarded one-electron Green function of the interacting system is given by the 
resolvent 

G(z) = [z-H cS (z)]- 1 (21) 
of an effective one-electron Hamiltonian with the matrix elements 

h r\,r'\'( z ) = #ra!r°a' + s a0) <5rr/ 5\\> • (22) 
The density of states (DOS) is defined as 

p{E) = -- Im Y, Gra,ra(£ + iO) . (23) 

71 A 

The set of equations has to be completed by an equation for the chemical potential fi 

r dEp(E) = n, (24) 

J — oo 

where n is the number of electrons per one site. 
The spectral density 

A(k,E) = --ImJ2 G xx (k,E + iO) (25) 

71 A 

is expressed in terms of the lattice Fourier transform of Grr'(z). 

NUMERICAL IMPLEMENTATION 

The formalism developed in the previous section was applied to the case of the 
paramagnetic fcc-Ni and the paramagnetic bcc-Fe which belong to weakly interacting 
systems. There is no exact way of determining the Hubbard interaction parameter U 
for the interaction among the d-electrons. One can find a broad variety of values in the 
literature. In the present calculations we have adopted U = 0.18 Ry for the paramagnetic 
fcc-Ni which emphasizes its atomic-like character, while the value U = 0.1 Ry used for the 
paramagnetic bcc-Fe is close to the value derived from experiment (see Ref. ||15|). The 
ratio U/w is larger for Ni as compared to Fe because the Ni bandwidth is smaller and 
the corresponding U is also larger. Therefore many-body effects will be more pronounced 
in Ni metal as compared to Fe. The smaller value of U/w in iron is, however, partly 
offset by a larger number of holes so that one can still expect non-negligible influence 
of many-body effects on one-particle states. The parameters of the one-particle part of 
the MBHH were obtained from the TB-LMTO assuming experimental lattice constants 
and the Ceper ley- Alder form of the exchange-correlation potential. We note that the TB- 
LMTO method employs the so-called atomic sphere approximation (ASA) which gives 
an accurate description of ground state properties of transition metals, their alloys and 
surfaces. For more details we refer to a recent book, Ref. |14] . 



The elements of the site-diagonal Green function corresponding to the one-particle 
part of the MBHH, Grr(-E), needed for the many-body calculations, are obtained from 
the site-diagonal elements of the auxiliary Green function, <7rr(z), z = E + iO, as 

G-r.\,rx'( z ) = [(z — H LMTO ) 1 ]ra,ra / 

= Xr\(z) 8x\> + ^kx(z) g-Rx,nx'(z) finx'(z) , (26) 
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where 

9nx,nAz) = ^J2[(V(z) - S(k))-V • (27) 
iV k 



In (p7[) , the sum runs over the Brillouin zone, N denotes the number of lattice sites, and 
the site-diagonal quantities Ara(z) and fin\(z) are functions of the potential parameters 



Cr, Ar, and 7r (see Ref. |L4| for more details). Finally, Vnx(z) is the potential function 



which is expressed in terms of potential parameters as 

Ara + (7ra - Prx){z - Cr A ) 

The calculations are performed first along a line in the complex energy plane (z = E +ie, 
|e| = 0.01 Ry), and the results are analytically deconvoluted back to the real axis fT8| . 



In this way, the poles present in Eq. ( P7j) are avoided. Once the imaginary part of 
Gra,ra(-E'+^0) is known, the selfenergy is determined from Eqs. (|T5|) - ClTTf) . This procedure 
is non-selfconsistent with respect to the Green function and the selfenergy. However, we 
determine the Fermi energy from (||) in a selfconsistent manner in order that the number of 
particles is conserved. In the last step we evaluate the quantities of interest, in particular, 
the density of states p(E), Eq. (|23|), and the spectral density A(k,E), Eq. (f25|). To this 
end we need to perform calculations indicated above but now for the Hamiltonian H eS , 
Eq. fl22"|) , rather then for ff hMTO _ i n other words, the potential function (p8|) is substituted 
by a function of the same form in which Cra is replaced by 

Cnx(z) = C RL + Sra(^) • (29) 

The calculations are again performed along the line z = E + it in the complex energy 
plane and then deconvoluted to the real axis. The selfenergy in the complex energy plane 
is obtained from its imaginary part on the real axis determined above using dispersion 
relation (||). 

We have used 2480 and 1360 k-points in the irreducible Brillouin zone to perform 



the integrations in Eq. (|27|) for the fee and bec cases, respectively. We note that such a 
high number of k-points is only needed to reduce the oscillations in the DOS tails which 
are due to the discretization of k-space, while for the determination of the selfenergy it 
is sufficient to use much smaller number of k-points. We have verified that the Fermi 
energy and the selfenergy obtained using 280 and 240 k-points give essentially the same 
result. The step in the energy was 0.005 Ry. A typical run for 750 energy points and 2500 
k-points requires about 5 minutes on a medium workstation and the selfconsistency with 
respect to the Fermi energy is achieved typically after 5 iterations in the weak-interaction 
case. 



RESULTS AND DISCUSSION 

In Fig. 2 we compare the DOS of the paramagnetic fcc-Ni evaluated in the non- 
selfconsistent SOPT and in the non-selfconsistent TMA. In each case the results are 
compared with the corresponding LDA-DOSs. The following points are to be mentioned: 
(i) we observe a narrowing of the one-particle (LDA) bands, which is in agreement with 
experimental data [[?j]; (ii) we note sharp structures close to the Fermi energy which are 



7 



only weakly influenced by electron interactions. The DOS at the Fermi energy is thus 
rather similar to that obtained in the LDA. Note that a large value of the DOS at the 
Fermi level is a precursor of a magnetic state at T = K (the Stoner criterion) which 
does not seem to be influenced by electron interactions; 




-0.8 -0.6 -0.4 -0.2 0.0 0.2 -0.8 -0.6 -0.4 -0.2 0.0 0.2 

E - ji (Ry) E - \i (Ry) 

Figure 2. Densites of states for a paramagnetic fcc-Ni within the LDA (dashed lines) and non- 
selfconsistent 2nd-order perturbation theory (left frame, full line) and non-sclfconsistent T-matrix ap- 
proximation (right frame, full line) for the Hubbard parameter U = 0.18 Ry. 

(iii) the 3d — 3d electron interactions in a partly filled band also give rise to a broad 
dispersionless satellite which appears below the Fermi energy. This satellite is positioned 
at about 6 eV below the Fermi energy in the SOPT in agreement with experiment M. It 
should be noted, however, that its position depends on the value of the interaction param- 
eter U used in the calculations. The results for the SOPT are in a reasonable agreement 
with similar calculations of Steiner et al. |fL5| ; (iv) the effect of electron interactions for 
the same value of the Hubbard parameter U is stronger in the TMA. In particular, the 
satellite structure is more pronounced and shifted deeper below the Fermi energy. This 
is due to the denominator of the T-matrix (see Eq. [16]) which is missing in the SOPT 
expression. 

We have performed similar calculations also for the RPA and the GWA. The results 
for the RPA are very similar to the results of the SOPT. The reason is that the fc-order 
diagram of the RPA is of order of 0(n£), where is the number of holes per site which 
is small for Ni so that the first term (k = 2) in the perturbation series, which is identical 
to the SOPT diagram, dominates. On the other hand, the effects of electron interactions 
(band narrowing and the shift of the satellite below the Fermi energy) are strongest within 
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the GWA. 

A deeper insight into the nature of the electron states influenced by electron-electron 
interactions could be obtained from spectral densities, or k-resolved densities of states. 
We note that the DOS and the spectral densities are, with exception of the transition 
matrix elements, proportional to the angle-integrated and angle-resolved photoemission 
spectra, respectively. Despite the fact that the effect of the transition matrix elements 
may be important, such quantities are still very helpful in understanding the influence of 
electron interactions on the electron states. In Fig. 3 we present the spectral densities 
of Ni metal calculated within the SOPT for k=(0,0,0) (D-point) and for k = ^(1,0,0) 
(X-point) in the BZ of the fcc-lattice together with the LDA spectral 



i i i i i i i i i r 
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Figure 3. Spectral densities for a paramagnetic fcc-Ni within the LDA (dashed lines) and non- 
selfconsistent 2nd-order perturbation theory (full lines) at k = T (left frame) and at k=X (right frame) 
for Hubbard parameter U = 0.18 Ry. 

densities. In calculating the spectral densities we added a small imaginary part (0.005 Ry) 
to the energy to avoid too sharp peaks, in particular for the LDA curves with no damping 
due to the selfenergy. Three LDA peaks for k = V correspond to three bands, one of 
the s-character at E ps —0.75 Ry and two other corresponding to the d-states around 
E —0.2 Ry. At the X-point we have in the LDA two peaks around E ps —0.3 Ry 
corresponding to the ci-states of which the lower one has a strong admixture of the s- 
states, the other two rf-bands at the Fermi energy (not resolved here because of their 
closeness and the finite imaginary part of the energy used in the calculations), and one 
sp-band above the Fermi energy (see Fig. 4). The influence of electron interactions results 
in a shift and broadening of the electron states, but the effect is anisotropic with respect 
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to k-vector and it depends also strongly on the energy region. 
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Figure 4. The LDA-bandstructures of paramagnetic fcc-Ni (left frame) and of paramagnetic bcc-Fe 
(right frame) along the lines of high symmetry in reciprocal space. 
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Figure 5. Real (full lines) and imaginary (dashed lines) parts of the selfenergy for a paramagnetic fcc- 
Ni for the non-selfconsistent 2nd-order perturbation theory: t2g- s y mme t r y (left frame) and e^-symmetry 
(right frame) for Hubbard parameter U = 0.18 Ry. 



To understand such a behavior in more detail, we plot in Fig. 5 the real and imaginary 
parts of the selfenergy H\(E). Note that our approximations obey the Luttinger theorem, 
i.e., ImEx(E) oc (E — fi) 2 in the neighborhood of the Fermi energy. Due to the cubic 
symmetry there are two different selfenergies, namely, T, t2g {E) and H e (E), corresponding 



10 



to the three- and two-dimensional representations for the d-states, respectively. Another 
point to be mentioned is that |ImE t2 (E)\ < |ImS e (E)\, because the e g band is narrower 
than the t<i g band, so that the ratio U/w is larger for e 9 -states and hence the influence of 
electron interactions is also stronger for e s -states. The states above the Fermi energy are 
only weakly damped by electron interactions and are shifted almost rigidly i.e., S Q (£') 
Re H a (E) ~ const., (a = t 2g , e g ). The same concerns the low-lying s-states at the bottom 
of the band (k = T) which are influenced by electron-electron interactions among the 
ci-electrons only indirectly, via a hybridization with the d-states. In both cases, the peaks 
have a Lorentzian shape indicating well-behaved quasiparticles. The same holds for the 
d-states close to the Fermi energy both for k = T and k =X. On the other hand, the 
d-states well below the Fermi energy, particularly in the energy region where |ImE Q ,(£')| 
have their maxima (E ~ —0.4 Ry), are strongly influenced by electron interactions. 
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Figure 6. Densites of states (left frame) and spectral densities at k=H (right frame) for a paramagnetic 
bcc-Fe within the LDA (dashed lines) and the non-selfconsistcnt 2nd-order perturbation theory (full lines) 
for Hubbard parameter U = 0.1 Ry. 



This is illustrated by spectral densities for k=X at energies around E m —0.3 Ry. The ex- 
tra spectral weight in the energy region E G (—0.6, —0.2) Ry which is missing within the 
LDA, gives rise to a satellite structure observed in the DOS. A strongly non-Lorentzian 
behavior of spectral densities in this energy region indicates a breakdown of the quasipar- 
ticle behavior. Note a particularly strong influence on the lower ci-band at k = X and a 
weaker influence on states at k = T, where the d-states are closer to the Fermi energy. We 
note that along the line T-X, the band-like features move with k-vector, while the extra 
spectral weight connected with the satellite structure, remains localized in the BZ. The 
change of the sign of Re£ a (-E) means that the states below and above E ~ —0.4 Ry are 
shifted in opposite directions. The band narrowing can be explained by a negative slope 
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of KeH a {E). The results for k =X agree reasonably well with the results of the SOPT 
calculations |T5| despite the fact that in the cited paper the authors studied a ferromag- 
netic Ni. Its exchange splitting is rather small. We have found that behavior of spectral 
densities in other approximations (TMA, RPA, and GWA) is qualitatively similar and 
thus need not be discussed separately here. 

We now briefly discuss the results for a bcc Fe to see differences due to a different 
structure, smaller U/w ratio, and a lower electron concentration. The results for the DOS 
and the spectral density at k= (1,0,0) corresponding to the H-point (see Fig. 4) are plotted 
in Fig. 6 for the SOPT. The results are in a reasonable agreement with those of paper |15| , 
although one has to keep in mind that our value of U is slightly larger and calculations are 
performed for a paramagnetic rather than for a ferromagnetic state. In the LDA, a bcc Fe 
exhibits well-known bonding and antibonding peaks with a high value of the DOS at the 
Fermi energy. This feature indicating a strong tendency of a paramagnetic Fe towards a 
magnetic instability (Stoner criterion) is preserved in the SOPT. On the other hand, the 
peaks below the Fermi energy are smoothed by many-body effects (finite imaginary part 
of the selfenergy) and due to a larger number of holes (partly filled <i-band) also states 
above the Fermi energy are now influenced. This is better seen in the spectral density 
evaluated at k=H. In particular, the LDA peak above the Fermi energy (see Fig. 6) is 
now much more broadened than similar peaks in the case of Ni. The other peak below the 
Fermi energy is influenced by electron interactions more strongly showing a weak satellite 
feature around E < —0.4 Ry. 



CONCLUSIONS AND OUTLOOK 

In summary, we have developed a scheme that allows to determine electronic prop- 
erties of correlated solids which could be described by a multiband Hubbard Hamiltonian 
and which belong to the weak- interaction case, U/w < 1. Our approach describes the one- 
electron part of the many-body Hamiltonian in terms of a TB-LMTO Hamiltonian. The 
TB-LMTO is the LDA-based technique which describes reasonably well the ground-state 
properties of solids. The present formalism is particularly suitable for future extensions 
to random alloys and their surfaces. This represents a great advantage of the present 
formulation in comparison to similar constructions found in the literature. The many- 
body part is described by a non-selfconsistent FLEX-type method, but the Fermi energy 
is determined selfconsistently. There are strong arguments [19| that such an approach 



is in fact more correct than simple selfconsistent treatments using dressed rather than 
bare propagators, because within an internally consistent theory one has to renormal- 
ize simultaneously propagators and vertices. We are working on such a generalization 
of single-channel approximations to include a properly selected subclass of parquet dia- 



grams p0[ . In the future, we plan to implement it into the present formalism. The basic 
approximation adopted in the present paper is a local character of the selfenergy. This 
limitation which is reasonably justified for transition metals is also motivated by future 
applications to alloys and their surfaces as it is consistent with the treatment of the sub- 
stitutional disorder using the CPA. On the other hand, the assumptions concerning the 
paramagnetic state, and, more importantly, the simplified form of the electron-electron 
interaction reduced here to a single value of the Hubbard parameter U is not essential 
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and can easily be lifted. We employ retarded quantities evaluated along a line in the com- 
plex energy plane coupled with an analytic deconvolution to the real axis and combined 
with the use of dispersion relations for the determination of the analytic function from 
its imaginary part. This is a powerful alternative to the conventional approach based on 
the use of the causal Green functions in which the energy integrals are replaced by finite 
sums over Matsubara energies. The knowlegde of the electron Green functions then allows 
one not only to evaluate one-electron properties, but also a future extension to transport 
properties of alloys and multilayers. 
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